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SUMMARY 

A comparative experimental investigation of the unstable operating 
modes of a centrifugal compressor is made. Hot wire measurements, in 
two different compressor configurations, allowed to distinguish between 
the different unsteady flow pattern. Impeller and/or diffuser rotating 
stall has been observed, depending on the flow conditions. The measu- 
red relative rotational speed of this perturbation is cross-checked 
with other experimental data and it is shown that the rotational speed 
is strongly dependent on the type of rotating stall . The diffuser 
absolute inlet flow angle at the onset of diffuser rotating stall 
agrees well with the value predicted by an existing stability criterion. 

INTRODUCTION 

Self excited vibrations of turbomachinery components, due to aero- 
dynamic forces, can be divided into two groups. A first one concerns 
the vibrations where the destabilizing forces are function of the 
whirling motion of the rotor. The initial displacement of the shaft 
influences the flow conditions in the labyrinth seals or creates 
an asymmetric blade loading on the impeller. The corresponding forces 
sustain the initial perturbation because of the negative damping. 

A second type concerns the vibrations where the main flow is de- 
stabilized by hydrodynamic or viscous forces which are independent of 
the resulting mechanical vibration. They are known as rotating stall 
and surge. The last ones have also an important influence on compres- 
sor performance. 

At low pressure levels, the corresponding forces are often Too 
small to have any significant influence on the shaft vibrations, but 
sometimes result in blade failures. However, at high pressure levels 
the excitation becomes much larger and very often limits the range of 
operation because of vibrational problems. 

This paper is limited to the study of the second type of desta- 
bilization in centrifugal compressors with vaneless diffusers. 
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Hot wire measurements have been performed in two different com- 
pressor configurations in order to investigate the different unsteady 
flow pattern. 

SYMBOLS 

b diffuser width 

Cp static pressure rise coefficient 

i number of stall cells 

P static pressure 

R radius 

Re[j Reynolds number, based on diffuser absolute inlet velocity 

and diffuser width 

U peripheral velocity 

V absolute velocity 

a absolute flow angle, measured from meridional plane 

B blade angle measured from meridional plane 

(|) mass flow coefficient 

p static density 

u slip factor 

tOg absolute angular velocity of stall cells 

n angular velocity of the impeller 

Subscripts 

0 settling chamber 

1 inducer inlet channel 

2 diffuser inlet 

3 diffuser outlet 

c critical 
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DESCRIPTION OF TEST FACILITY AND INSTRUMENTATION 

The experiments have been performed in the open loop test facility 
(Fiq. la) of the State University of Gent. Mass flow is fed into the 
Uttling chamber by a normalized tube which allows for 
and mass flow measurement. Inlet guide vanes can be 

suction pipe if prerotation is required. Flow is discharged at atmos- 
pheric pressure after a 180° turning at the diffuser exit. 

All measurements have been performed with a vaneless diffuser. 
Geometrical parameters of the unshrouded impeller (Fig. 2) are listed 
in Table 1 . 

Besides the classical pressure and temperature probes, required 
for overall performance measurements, four hot wire probes are ^ 
installed. The hot wires are mounted perpendicular ^ 

direction. One is installed at impeller (II), 

inlet (ID) and two at R/Rz = 1-3 inside the diffuser (MDl and MD2) at 

42° circumferential distance (Fig. lb). The 1 channel 
visualized on an oscilloscope and processed with a HP 3582A two channel 
FouHeJ spectrum analyser. This allows to define the power spectrum, 
coh6rGnc6 function dnd phdS6 trdnsfor function- 

Perturbations rotating in the diffuser will result in 
signals at MDl and MD2 with period tz . Because of the 
in circumferential position there will be a time lag ti which allows 
to calculate the number of stall cells i, and the rotational speed 

i = 360/42 * ti/tz (i must be integer) 


0 ), 


= 2tt/ ( i T 2 ) • 


TEST RESULTS WITHOUT PREROTATION VANES 

The first series of measurements have been obtained with an axial 
flow at compressor inlet. The overall 

figure 3. The static pressure rise is not related to the mass flow 
coefficient 4.9 but to the diffuser inlet flow angle aa • The last 
parameter is used here because of its important influence on diffuser 
stability and is defined by 


tg a 2 ” u/4^2* 

The slip factor is calculated from correlations and is checked against 
the measured impeller temperature rise. Increasing values 0 02 

corrGspond to decrGSSing iridss flow. 

The static pressure rise is non dimensional ized by P 0 U 2/2 to 
become independent of RPM. 


All curves show a maximum at 02 = 76° which coincides with the on- 
set of unstable flow. At higher RPM, these curves do not continue to 
the left because of surge. At low RPM, surge occurs only at 02 
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which leaves a large range of flow conditions where unsteady flow can 
be observed. The hot wire measurements are therefore made at 2000 RPM 
for aj values between 76° and 87°. 

• 1 traces shown on figure 4a are measured at rotor 

inlet (II) and diffuser inlet (ID) for aa = 76°. The relative ampli- 
tude of the oscillations, defined by ( V^ax'^'mi n ’ ^'^O'^nts to .14 

at rotor inlet and .23 at diffuser inlet. The power spectrum and co- 
herence function, shown on figures 4b and 4c, clearly indicate an 
instability at 27.2 Hz. Inside the diffuser, at R/R 2 = 1.3, the flow 
IS steady as shown on the power spectrum of probe MDl (Fig 4d). One 
can therefore conclude that these instabilities are due to impeller 

stall, which can be explained by the high positive incidence 
(12 ) at inducer tip for this operating point. 

Measurements at a 2 = 78° (Fig. 5a) do not only show an increase 
in amplitude to .37 at impeller inlet (II), but also an extension of 
pese perturbations downstream inside the diffuser (MD). The hot wire 
traces MDl and MD2 (Fig. 5b) have a relative amplitude of .13. The 
phase lag between both signals indicates the presence on one stall 
cell rotating at 11% of the rotor speed. The power spectrum and co- 
herence function of the II and MD signal (Figs. 5c, 5d) reveal a 
decrease of the basic frequency to 25.6 Hz and the presence of impor- 
tant harmonics at 52 and 77 Hz. 

Figures 6a and 6b are measured at a 2 = 81° and indicate a further 
increase in amplitude to 200% at rotor inlet (II) and 18% inside the 
diffuser (MD). The basic frequency is now 23.2 Hz and wc/n = .70 The 
harmonics have almost disappeared at rotor inlet but are unchanged 
inside the diffuser. 

Low frequency oscillations completely suppress impeller rotating 
«2 = 87 . The hot wire traces of figure 7a seem to indicate 
that the frequency at rotor inlet (II) is twice the frequency of 
oscillations inside the diffuser. However, hot wires read only the 

perpendicular to the wire and are insensitive to a 
180 change in flow direction. The real flow velocity at rotor inlet 
changes sign (return flow) and is indicated by the dashed line. Both 
signals are in phase at 6.4 Hz and indicate compressor surge The 
power spectrum on figure 7b shows the presence of several harmonic 
frequencies as well at rotor inlet as inside the diffuser. 

This series of measurements show only one type of rotating stall 

Z'ilniill oscillations are^radually inc^ea^•ng 

in amplitude and extending downstream inside the diffuser. They are 
due to one stall cell rotating at high relative velocity. 

TEST RESULTS WITH INLET PREROTATION VANES 

/n f According to a stability criterion, developed by the authors 
(Ref. 1), vaneless diffuser rotating stall should occur at a, = 77° 
However, only impeller rotating stall has been observed in previous' 
experiments. A possible explanation for this disagreement could be 
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that diffuser rotating stall is completely suppressed by the impeller 
rotating stall if the latter one starts first. 

In order to verify this hypothesis, inlet guide vanes have been 
installed in the suction pipe, in order to create 60° 

the impeller inlet. Preliminary calculations indicated that this^ 

would delay impeller rotating stall at least to «2 - . In this way 

a sufficiently broad margin for diffuser rotating stall would be 

created . 

The overall compressor performance map is ^hown on figure 8 
Comoared to previous results, one can observe a 6 shift to the lett 

0 ? im^ener characteristics. The maximum ^P^nn'^'T^dm s 
has slightly decreased because of the prerotation. The diffuser 

pr6ssur6 rise is almost unchanged. 

First instabilities are observed in the diffuser at U 2 = 77° 
rPia 9al The frequency is 17.6 Hz (Fig. 9b) with a coherence of 
5^^^(Figl ’9c) . Frol the 132° phase shift between MDl and MD2, one 
can calculate that three stall cells are rotating at . 

These instabilities are not observed at the inducer inlet (Fig. 9d) 

conclude that they are due to diffuser rotating 

stal 1 . 

These instabilities have a maximum amplitude of .11 at a 2 .= 78° 
but completely disappear at «2 = 80°. They suddenly reappear in the 
diffuser at «2 = 80°5 with a relative amplitude of .14 and a cohe 
rence of 98%. 

With the achieved prerotation, the critical flow angl e for impel - 
ler rotating stall is predicted at «2 = 82°. This is experimentally 
verified by the power spectra obtained at az T o. The power peak 
at 23.2 Hz is measured at rotor inlet and inside diffuse . _ 

phase angle of 52° indicates that one stall ce 1 rotates at wg/n - 
^673. These results are very similar to "^all 

observed with zero prerotation. The power peak at 19.2 Hz observed 
only inside the diffuser and the phase angle between MDl 
cates that three stall cells are rotating at 

to the measurements at lower az values, one can conclude that this 
diffuser rotating stall. 

When closing the throttle valve further to az = J ^ 
rotating stall remains unchanged, but impeller rotating stall 
(23.2 Hz) disappears completely (Fig. H). 

A new type of unstable flow, with a frequency of 8.8 Hz is 
observed at U 2 = 83»6 (Fig. 12a). These instabilities are observed 

only in the diffuser and correspond to a °%J’*^*^!fheL''high 

ctail Two Stall cells are rotating at ~ 13.2%. Ihese nign 

and^low frequency diffuser rotating stall patterns do not occur simul- 
taneously b3t in a rather intermittent mode. This relative amplitude 
increases further to .22 at a 2 - 84 7 (Fig. 12b). 
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The low frequency oscillations (6.4 Hz), observed at a, = 86° 
beared" a1'lhirpo?n?.‘° ""toting stall has d1sap- 

COMPARISON OF THE TWO TEST SERIES 

^ observed at impeller inlet and 

series, the onset of this type of 
stall coincides with a 12 positive incidence at impeller tip. Pre- 

laH nSw C0e???c?enu°"' instability to lower 

start?^It'^a^''-'"77o^^'^Th \ Observed in the diffuser and 

next chapter ~ value will be compared to other data in the 

Both types of rotating stall can exist simultaneously if diffuser 
i?^nf ^ occurs first. This absolute frequency and amplitude 
IS of the same order of magnitude and measurements with one sinale 
probe do not allow to distinguish between them. However, the re ative 
rotational speed is quite different as shown in figure iL 

7 rotating stall has a relative rotational speed of about 

''® ocity seems to be unaffected by prerotation. The rela- 
^in tn of diffuser rotating stall is of the order of 

. u to only. This difference is in agreement with experimental 

^ Howard (Ref. 2). Misuki et al.Vef. 3) 
and Abdelhamid (Ref . 4) . This difference could therefore be used to 
distinguish between impeller and diffuser rotating stall. However 
the definition of relative rotational speed not only requires t^ 

measure the absolute frequency but also the phase difference between 
two circumferential positions. ueuween 


THEORETICAL PREDICTIONS OF VANELESS DIFFUSER 
ROTATING STALL 

Several stabil i ty criteria for vaneless diffusers have been 
developed and published in the literature. They are based on dif- 
ferent assumptions and result in contradictory conclusions. 

/D stability criterion for vaneless diffusers of C. Rodgers 

(Ref. 5) IS based on theoretical considerations concerning variation 
of Cp in function of ^2 • The critical value of azc, separating the 
zone of stable operation from the unstable zone, is function Sf the 
non dimensional ized diffuser width b2/R2 and radius ratio R3/R2. 

/D ^ conclusion is similar to the one obtained by W. Jansen 

(Ket. 6) His criterion is based on detailed calculations of the 
VISCOUS flow in vaneless diffusers. Besides diffuser width (b?/Ro) 

and flow angle («2 ) , stability is now also^ a func 
tion of the Reynolds number. Both theories suggest that the critical 
flow angle a2c is decreasing for decreasing diffuser width. This 
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would mean that for a given impeller, an unstable vaneless diffuser 
cannot always be stabilized by decreasing the diffuser width. 


A similar but more precise calculation by Senoo (Ref. 7), however, 
results in a criterion which shows an opposite trend in terms of 02 c 
as a function of diffuser width. The critical f ^ 

increasing for decreasing diffuser width. An unstable diffuser can 
therefore be stabilized by decreasing diffuser width. 

also allows to draw important conclusions in terms of diffuser radius 
ratio R 3 /R 2 3"^^ about the influence of non uniform flow at diffuser 

inlet. 


The stability criterion of A. Abdelhamid (Ref. 8) is based on 
rotor-diffuser interaction and does not account for viscous effects. 
This criterion relates stall operation to the pressure and tangential 
velocity fluctuations non dimensional ized by radial velocity fluctua- 
tions. This criterion is not easily applicable and stability is inde 
pendent of diffuser width. 


In order to eliminate the contradiction between the different 
criteria, the authors (Ref. 1) have performed some unsteady flow 
measurements on different impeller-diffuser combinations. These 
results and other experimental data from the 1 1 terature al 1 owed to 
derive the following correlation. All values of asc within a 

relatively narrow band as shown on figure 14 after they have been 
corrected for Reynolds number influence according to figure 15. 
Points 1 to 5 are experimental data obtained by the authors. Points 
6 to 16 are derived from published data. 


These experimental data indicate an increase of a 3 c for decreasing 
diffuser width and agree very well with the predictions of Senoo. His 
theoretical results do not indicate an influence of Reynolds number on 
diffuser stability for values of b 2 /R 2 < .10. However, the experimen- 
tally observed variation of a.2c function of Reynolds number is . 

agreement with the predictions of Jansen (Acx 2 = 6 = between Reb^ - l 

and Rob^ = 10^) . 

The flow angles 1 to 16, are calculated from measured temperature 
rise and slip factor correlations. Therefore, some doubt could exist 
about the precision by which these angles are defined. Somespecial 
tests have been performed in order to evaluate this error. 
made for one impeller-diffuser combination at two different Reynolds 
numbers. Flow angles at diffuser inlet are measured by a crossed hot 
wire. The experimental values shown by points 17 and 18 are already 
corrected for Reynolds number influence. 


The observed differences with the theoretical curve are within an 
acceptable tolerance and therefore provide an additional verification 
of this correlation. 


Point 19 is the value obtained from the present test series and 
confirms the stability limit for diffuser rotating stall. 
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CONCLUSION 


JIJOW lo t:^‘;j5ufsh^between 

speld irab^r 70 ^"* incidence and the relative rotJtlonll 

inlet air angle’and Reyno?d1^™ber\^%1’ra’re^H;r;St%?(onirsp^ 

nhan experiments provide a criterion to distinguish between both 

aHoufd1r::3if?:S 

Flow conditions at the onset of diffuser rotatinn ctaii h=.wrv k 
compared with a criterion for diffuser stabfutTaJd agree ie^y wen"" 
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TABLE 1. - GEOMETRICAL CHARACTERISTICS OF THE COMPRESSOR 
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Figure 13. - Relative rotational speed. 
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Figure 14. - Critical <nlet flow angle complement versus nond imen s i on a 1 
diffuser width. 



- critical inlet flow angle correction factor versus Reynolds 
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